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Quantum confinement is known to infiuence fermionic condensates, resulting 
in quantum-size oscillations of superfiuid/superconducting properties. Here we 
show that the impact of quantum-size eff"ects is even more dramatic. Under re- 
alistic conditions, a significant phase-space reconfiguration induced by quantum- 
size eff"ects opens a quasi-molecule channel in the fermionic pairing so that the 
condensed pairs exhibit features typical of a molecular state. As an illustration 
we consider a quasi-one-dimensional fermionic condensate, as realized, e.g., in 
cigar-shaped atomic Fermi gases or superconducting quantum wires. In this case 
the transverse quantization of the particle motion favors pairing through a co- 
herent superposition of quantum channels that are formed due to the grouping 
of single-particle levels into a series of well distinguished subbands. Whenever 
the bottom of a subband approaches the Fermi level, the longitudinal spatial dis- 
tribution of fermions in a condensed pair becomes strongly localized within the 
corresponding quantum channel. The fermionic pairs in this channel resemble 
molecules with bosonic character. 

In 1960's Blatt and Thompsoi>i'^ calculated the energy gap of an ultrathin superconduct- 
ing slab as function of its thickness and found a series of pronounced peaks. They called 
these peaks shape superconducting resonances. The physics behind this result is usually 
understood as follows. The spectrum describing the electron motion in the direction per- 
pendicular to the film is quantized whereas a quasi-free electron motion is assumed in the 
direction parallel to the film. Due to the transverse quantization, the conduction band splits 
up into a series of single-particle subbands. The lower edges (bottoms) of such subbands 
are determined by the perpendicular discrete electron levels and, so, move in energy with 
changing thickness d (scaling as ~ l/rf^). Unlike the subbands, the Fermi level exhibits a 
significantly less pronounced size-dependent shift. So, the bottoms of the subbands pass 
through the Fermi surface one by one when increasing the thickness. Each time when the 
bottom of a subband approaches the Fermi surface, the density of single-particle states in 
the vicinity of the Fermi level increases, resulting in an enhancement of basic supercon- 
ducting quantities, e.g., the critical temperature, the order parameter and the excitation 
gap. Thus, the formation of single-particle subbands due to quantum confinement mani- 
fests itself through quantum-size oscillations of the superconducting properties driven by a 
series of shape superconducting resonances (below they are also referred to as quantum-size 
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superconducting/superfluid resonances or size- dependent superconducting/superfluid reso- 
nances) . 

For several decades after the paper by Blatt and Thompson, only few experimental groups 
reported possible signatures of quantum-size oscillations in superconducting films'^"-. How- 
ever, poor sample purity and significantly limited size control prevented definite conclusions. 
Blatt and Thompson also argued (see Refs. flfl) that similar physics could be expected for 
nucleon pairing in atomic nuclei as a consequence of the shell structure of the single-particle 
spectrum. This expectation was found to be in agreement with experimental data^ and, for 
the next 40 years, atomic nuclei were the only system where quantum-size effects on the 
BCS fermionic pairing were investigated both experimentally and theoretically (for more 
details, see Ref.0 and references therein). Advances in microstructuring of superconduc- 
tors in the 90 's renewed interest in quantum-size effects in superconducting systems. The 
work of Blatt and Thompson was extended by Bianconi, Perali and coworkers^iii^ to quan- 
tum striped superconductors, where a sizeable amplification of the critical temperature and 
excitation gap was found in the presence of a shape superconducting resonance. Recent 
developments in nanofabrication have resulted in high-quality superconducting nanosys- 
tems, and the quantum-size oscillations of the critical superconducting temperature were 
eventually observed in single-crystalline atomically uniform in thickness Pb nanofilmsi^"— . 
Furthermore, the results of a numerical self-consistent solution of the Bogoliubov-de Gennes 
(BdG) equations in Ref. UM showed that the quantum-size superconducting resonances were 
responsible for a systematic thickness-dependent shift-up of Tc found in high-quality alu- 
minum and tin superconducting nanowires^^— . Very recently quantum-size oscillations of 
the excitation gap were also reported for tin superconducting nanograins^i. 

Cooling of trapped fermionic atoms down to ultra-low temperatures such that the atomic 
gas is Fermi degenerate (see, e.g., Ref. |22| ). resulted in another system promising for the 
investigation of quantum-size effects in fermionic condensates. In the case of trapped atomic 
Fermi gases the single-particle levels can be tuned by laser light or magnetic fields through, 
e.g., changing the spatial dimensions of the trap. In particular, for a pancake- or cigar-shaped 
geometry the perpendicular trapping frequency Uj_ is much larger than the parallel one a;||. 
In this case single-particle states form well distinguished subbands and the physical picture 
turns out to be similar to that in superconducting metallic nanowires and nanofilms, i.e., 
a quantum-size superfluid resonance can be expected when the bottom of a single-particle 
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subband is located in the vicinity of the Fermi level. In particular, quantum-size oscillations 
of the critical temperature and order parameter were recently calculated for a superfiuid 
Fermi gas confined in a quasi- two- dimensional (quasi-2D) harmonic trap^^. We note that 
fermionic condensates are now attainable even in the tight-confinement regime with only 
one perpendicular single-particle level below the Fermi energy^. Recently, experimental 
results on the first indications on quantum-size oscillations of the cloud-size aspect ratio 
in a harmonically trapped quasi-2D Fermi gas became available^^. Maintaining the atomic 
cloud at the lowest attainable temperature and at a low number of atoms, the group of 
Vale^ was able to detect the filling of the individual perpendicular levels, following in 
a quantitative way the dimensional crossover from 2D to 3D in a ^Li Fermi gas, while 
keeping the 3D character of the two-particle scattering. This work clearly demonstrates that 
the experimental conditions to explore such a dimensional crossover and the accompanying 
quantum-size effects are presently achievable with ultracold atoms. 

In this paper, we consider a quasi- ID fermionic condensate and show that the impact 
of the formation of single-particle subbands is more dramatic than simply resulting in os- 
cillations of Tc or the excitation energy gap. Significant reconfiguration of the phase space 
due to grouping of single-particle levels into a series of subbands (here quasi-lD means that 
we deal with more than one ID subband of single-particle states) results in opening of a 
quasi-molecule channel in the fermionic pairing each time when the bottom of a single- 
particle subband approaches the Fermi surface. Due to a "depletion" of the longitudinal 
Fermi motion in such a subband the longitudinal spatial distribution of fermions inside a 
condensed fermionic pair squeezes and becomes strongly localized. In this case the system 
behaves similar to the well-known crossover from the Bardeen-Cooper-Schrieffer (BCS) pairs 
to the Bose-Einstein condensation (BEC) of quasi-molecules driven through the Feshbach 
resonance in ultracold Fermi gases (see, e.g., Ref. |22|). However, very different from previous 
case, the effect of interest is now controlled by the quantization of the perpendicular par- 
ticle motion. By changing the transverse size of the system, one can change the energetic 
positions of the transverse single-particle discrete levels (bottoms of the single-particle sub- 
bands) with respect to the chemical potential fi. This leads to a significant redistribution of 
the kinetic energy between the transverse and longitudinal degrees of freedom in a subband 
and, so, to the above mentioned "depletion" of the longitudinal Fermi motion when the 
bottom of this subband approaches /i. Below, for illustration, we investigate a cigar-shaped 
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atomic Fermi gas and a superconducting quantum wire. Similar results are also expected 
for quasi-2D fermionic condensates. 



I. HARMONICALLY TRAPPED QUASI-ID SUPERFLUID FERMI GAS 

We begin with the BCS self-consistent equation for the spin-singlet s-wave pairing gap 
Ak of a 3D Fermi gas with an attractive interaction, i.e., 

1 1 



Ak = - ^ ^kk' Ak' tanh(^£;k'/2) 



(1) 



where Tk = fi?k'^ /2M is the single-particle dispersion, i?k = a/ (^k — /^)^ + ^k stands for the 
quasiparticle spectrum, /i is the chemical potential and /3 is the inverse temperature (/3 = 
1/kBT). The standard scattering length regularization is introduced in the gap equation 
(1) to avoid ultraviolet divergency in 3D. The interaction matrix element in Eq. ([1]) can be 
written in the form 

V^^> = -g [ dVl^kWn^k'WP, (2) 



where g = Anh^lal/M, with a < the 3D s-wave scattering length, and v'k(i') is the single- 
particle wave function that is proportional to e*^"" (plane waves) in 3D. 

When switching to a Fermi gas confined in a trap, the particle momentum label in Eqs. 
and ([2]) should be replaced by an index for discrete energy levels because the relevant wave 
functions are not plane waves but solutions of the one-particle Schrodinger equation for 
the corresponding confined geometry. In this section we consider harmonically trapped 
fermions in an axially symmetric confining potential M{ujj_p'^ + w^2;^)/2 (with u± ^ u)\\), 
where cylindrical coordinates are invoked. Due to the axial symmetry the relevant quantum 
numbers can be chosen as n = 0, 1, . . ., the radial quantum number; m = 0, ±1, ±2, . . ., 
the azimuthal quantum number; and j = 0, 1, . . ., the quantum number associated with the 
parallel (longitudinal) motion. For the aspect ratio of the trap we have l\\/l± ^ 1 (with 



/|l = ^JftjMiiJw and = ^yh/Muj_) and, so, using the 3D pseudopotential g = Anfi?\a\/ M 
seems problematic due to the possible ID character of the two-particle scattering. However, 
the aspect ratio is not actually the quantity that controls the character of scattering in 
a cigar-shaped many-particle system. The analysis of the two-particle scattering in an 



axially confined geometry, see, e.g., Ref. 1271, shows that the scattering rapidly becomes 



three dimensional when the number of contributing subbands (perpendicular levels) is more 



than 1. For parameters considered in this work this number varies from 6 to > 20. Thus 
we are far beyond the regime of the effectively ID scattering and, so, any issue related 
to the confinement-induced Feshbach resonance in a ID wave guide and the formation of 

n nn 

confinement-induced molecules in a ID Fermi gas (see, e.g., Refs. '22' andl28H31l) are beyond 
the scope of the present work. To avoid any confusion, we would like to note once again 
that quasi-lD in the present work means that more than one transverse single-particle level 
contribute to the basic physical quantities. 

By numerically solving Eqs. ([1]) (with k ^ u = {n,m,j}), we self-consistently calculate 
for a given temperature a set of pairing gaps A^. Then, the critical temperature can be 
found as the temperature above which only the trivial (A,^ = 0) solution to the problem 
exists. To obtain more information, e.g., concerning the spatial distribution of the pair 
condensate and the fermionic pairing correlations, we should take into account that Eq. ([1]) 
follows from the Bogoliubov-de Gennes equations provided that the particle- like and hole-like 



wave functions Uu{r) and Vu{r) are approximated as (see, e.g., Refs. I32lj33l) 



M^(r) = W^V9j,(r), v^{r) = V^ip^{r), (3) 
with (Uu and Vu are taken real) 

and V5i.(r) = {}„m{p,f)Xj{z), where 'dnm{PiV) ^'^^ Xj{^) are the eigenfunctions of the 2D 
(isotropic) and ID harmonic oscillators, respectively. This approximation accounts for the 
pairing of time reversed states, and, as known since the pioneering paper by Anderson^i^ such 
a simplification yields quite reasonable results in the presence of time reversal symmetry. 
Then, based on the Bogoliubov transformation of the field operators V't('^) ^^"^ ''^ii^) the 
quasiparticle quantum amplitudes 'yu,t and 7^,4., i.e., 

one can study the superfluid pair correlations and the spatial distribution of the pair con- 
densate. In particular, introducing the Cooper-pair "wave function" \E'(r,r') = (■?/'|(r)-?/'4^(r)), 
we find 

^(r, r') = Y ^nn.ij. r ), ^„„^(r, r') = 79„„(p, ^) Cm(p', i^nm{z, z'), (4) 



with 



^nmj tanh (^/3Enrnj 



(5) 



where Tnmj = ^±(1 + 2n + \m\) + hw\\{j + 1/2). It is worth noting that the regularization 
term l/{2Tnmj) appears in the parenthesis of Eq. (|5]) as a simple extension to spatially 
nonuniform systems of the regularization used in Eq. ([1]). It is important to note that we 
have verified that different possible regularization schemes lead to insignificant differences 
in the physical quantities, preserving the basic conclusions of the present work, see a more 
detailed discussion in the end of this section. The single-particle states with the same n and 
m form a quasi-continuum. Therefore, it is natural to introduce single-particle subbands 
denoted by (ra, m) (see Fig. 1 in Supplementary Information) and the subband dependent 
fermionic-pair "wave function" defined by Eq. We remark that treating \E'(r,r') as the 
wave function of a condensed pair of fermions goes back to the classical papers by Gor'kov^- 
and Bogoliubov^^ and is directly related to the conventional interpretation of the order 
parameter A(r) = (7^1/ (r, r) as the center-of-mass Cooper-pair wave function. However, there 
exist also other ways to introduce the fermionic-pair wave function (see, e.g., Ref. [371 and 



38l ). Nevertheless, all these variants were shown^ to result in the bulk Cooper-pair radius 



being proportional to the ratio of the Fermi velocity to the energy gap. In the present paper 
we are interested in the longitudinal Cooper-pair size in a quasi- ID fermionic condensate 
which can be calculated as 



{^{z-z'f\^) 



1/2 



(6) 



with \I' given by Eq. @j. Below we deal also with the subband-dependent longitudinal size of 
condensed fermionic pairs ^g""^^ which is controlled by ipnm{z-, z') and calculated from Eq. (j6]) 

with ^(r,r') replaced by lljnmiz, z'), i.e., ^^"""^ = [{ll)nm\{z - ^')^l^nm)/(^nm|V'nm)]^''^- 

For a numerical solution of Eq. ([T]) (with !<;—)■// = {n, m,j}), we consider a quasi-lD 
harmonically trapped mixture of ^Li atoms in the two lowest spin states, \F^mp) = 1 1/2, 1/2) 
and |l/2, —1/2). The pair interaction between these states can be significantly modified by 
an external magnetic field by means of the formation of a broad Feshbach resonance^S.. In 
particular, typical values of the s-wave scattering length a at the BCS side of this Feshbach 
resonance varies^ from —250 to — lOOnm dependent on the magnetic field B (theoretically 
a goes to — oo when approaching the point of the Feshbach resonance B = 0.83 kC). In 
Fig. [Dour numerical results are shown for a = —140, —180 and — 210nm [panels (a,b). 
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(c,d) and (e,f), respectively]. The chemical potential and the longitudinal frequency are 
kept constant, i.e., fi = h ■ 24kHz and hw\\ = 0.01/i {uj\\/2'k = 240 Hz), while the trapping 
frequency in the perpendicular direction varies in such a way that discrete single-particle 
levels for the perpendicular single-particle motion pass through the Fermi surface when the 
ratio fi/hu)j_ reaches 2, 3, 4 etc. In our cigar-shaped confining geometry uj\\ is taken much 
smaller than u± ~ fi/h. The particular value of u\\ in this regime appears to be of no 
importance (changes less than a few percent were found for u\\ = 2tt ■ 300Hz and 27r ■ 480Hz). 
The chemical potential is chosen such that experimentally accessible values of the particle 
density and the ratio T^/Tp are obtained, with Tp denoting the Fermi temperature (for more 
details, see the discussion at the end of this section). 

Figured] (a) shows our numerical results for Tc (in units of ft/ks) as a function of fi/hu± 
calculated for a = —140 nm. As seen, each time when a discrete transverse level crosses fi, 
is enhanced, i.e., a superfiuid quantum-size resonance occurs. As a result, Tc oscillates with 
changing w^. These oscillations are damped: their amplitude is reduced when uj±_ decreases 
{l± increases). Such a decay is typical of quantum-size oscillations of the basic physical 
quantities in quasi-lD and quasi-2D fermionic condensates^^— The reason is that the 
number of single-particle subbands making a contribution increases with l^. In particular, 
for fi/hLo± = 2 there are six subbands that are responsible for the formation of 99% of the 
condensate [subbands (n,m) = (0,0), (0,±1), (0, ±2) and (1,0), major contribution is due 
to {n,m) = (0,±1)]. For fi/hw± = 4 there are already 21 contributing subbands. The 
larger this number is, the less pronounced is the increase in the relevant density of states 
when a new perpendicular level crosses /i. Then, quantum-size resonances are weakened and 
the corresponding oscillations are finally washed out. Another important feature of such 
oscillations is that their amplitude is also reduced when the pair interaction is enhanced 
(i.e., \a\ increases). A stronger pair interaction results in larger pairing gaps. In this case 
contributions of the levels even far below or above /i become of importance. As a result, 
passages of perpendicular levels through the Fermi surface have less pronounced and rather 
smoothed effects. This is illustrated by Figs. [T] (c) and (e), where /c^Tc/yU is given versus 
fi/hujj_ for a = — ISOnm and a = — 210nm. Note that in the present work is evaluated 
at the mean-field level. On the other hand, it is well known that the mean field critical 
temperature through the BCS-BEC crossover (driven by a change in the coupling kpa) cor- 
responds to a characteristic temperature of the pair formation and pseudogap opening, while 
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the critical temperature for the superfluid transition is renormahzed by pair fluctuations. 
In the simplest approach, the renormahzed critical temperature can^ evaluated by using 
a non self- consistent t-matrix approach, as shown in Refs. 



23 



and 



24. For our choice of 



parameters we have —1.4 < {kpa)^^ < —0.9 (with kp the Fermi wavevector in the center of 
the trap, see the discussion at the end of this section). These values of kpa are between the 
BCS and the crossover regimes, where the suppression of the critical temperature due to the 
pair fluctuations is not larger than 16%, thanks also to the fact that in trapped systems the 
fluctuation region shrinks with respect to the homogeneous case. Therefore, the amplitude 
and period of the oscillations of induced by quantum conflnement as shown in Fig. [H will 
not be importantly modifled by fluctuations. 

Based on the above results for T^, one can also anticipate similar oscillations of the longi- 
tudinal BCS coherence length ,^o- These oscillations are shown (red solid circles) in Figs.[^b), 
(d) and (f) for a = —140, —180, and — 210nm, respectively. Naive expectation based on the 
ordinary BCS picture suggests that oc hvp/{kBTc), with the Fermi velocity Vp = kkp/rrie. 
However, this expectation does not match our numerical results. As seen from Fig. [T](a), the 
coherence length drops by a factor of 2.7 when fi/huj± increases from 1.8 to 2.1. While 
according to the above naive estimation, is expected to fall only by about 40% (the critical 
temperature increases by about 50% with almost negligible change of vp). Similar difference 
can be found between panels (d) and (f). The reason for this difference is the formation of 
the multi-subband structure. For illustrative purposes, the subband-dependent coherence 
length ^o"™^ is also shown in Figs. HJ^b), (d) and (f) for (n, m) = (0,0), (0, ±1), (0, ±2) and 
(0, ±3). At /i = 2hw± the bottoms of the two degenerate subbands with {n,m) = (0, ±1) 
pass through the Fermi surface and the corresponding quantum-size resonance develops (a 
subband whose bottom is in the vicinity of the Fermi surface is referred to as resonant). At 
this point subbands (0, ±1) make the major contribution to the basic physical quantities. 
In particular, their common contribution to the integral of |^'(ri,r2)P is about 70 - 75%. 
In other words, we may state that almost 70 -75% of the fermionic pairs come from these 
subbands. As a result, being close to ^o°'°^ at /i < 1.8 hhj±, C,o changes its trend abruptly and 
approaches ^^'^^ = when the perpendicular level 2huj± crosses jj. For a = — 140nm, 

panel (b), this change is more dramatic because of a more pronounced difference between 
ande^^'^ When \a\ increases, ^o°'°^ and ^l)''^^^ become closer to one another and, as 
a result, the effect is weakened. A similar weakening also occurs for upper- level resonances. 
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i.e., associated with the perpendicular levels 3Huj±, Ahuj±, etc. This is mostly because the 
inter-subband energy spacing is reduced {hw^ decreases) so that the difference between the 
subband-dependent lengths becomes less and less pronounced. It is worth noting that the 
above naive estimation of the longitudinal Cooper-pair size, i.e., oc hvp/ikBTc), is not to- 
tally irrelevant: it yields reasonable results for ^q"™'' ^ subband with the bottom far below 
fi. For instance, from Fig. Wijo) one finds that ^0°'°'' reduces by about 30% when ^/huj± 
increases from 1.8 to 2.1, which is close to the drop of 40% following from the above sim- 
plified estimation. In addition, one can compare ^0°'°'* jj,/huj± = 1.8 (below the resonance 
point fi/hw± = 2.0) for different panels (different interaction strength) in Fig. [TJ As seen, 
this quantity scales approximately as the inverse critical temperature, which agrees with the 
estimate hvp/ikBTc). Now, let us check what happens with a subband whose bottom is in 
the vicinity of fi. Here the naive estimation of the longitudinal Cooper-pair size is not longer 
in agreement with our numerical results. In particular, ^^'^^^ taken at /i = 2huj± does not 
scale as 1/Tc when changing a (the bottoms of the subbands with {n,m) = (0, ±1) crosses 
the Fermi level at /i = 2fvjj±). Here a reasonable agreement can be achieved when assuming 
a less sensitive scaling, e.g., 1/a/7^. The same occurs for ^0°'^^'' and ^o^'"*^^^ at /i = 3huj± 
and Ahw±, respectively. However, scaling 1/a/T^ does not match at all when the bottom of 
a subband goes above the Fermi surface. As seen from panels (b), (d) and (f), ^nm's for 
[n, m) = (0, ±1), (0, ±2) and (0, ±3) are close to one another for fi < 2hjjj± and practically 
do not change with a. It is worth noting that ^0^'°'' and ^0"'^^'' exactly the same in 

spite of the fact that subbands (0, ±2) and (1,0) are degenerate. This difference is within 
several percent and appears due to a difference in the relevant interaction matrix elements, 
which results in ipi^o{p,(f) ^ '4'o,±2{p,f)- A similar difference appears between ^0°'^^^ and 
so 

To go in more detail about the subband-dependent fermionic pairing, we consider how 
tljnm{z,z') decays with increasing \z — z'\ (the characteristic length for this decay is ^0"'"'') 
for different energetic positions of the bottom of the corresponding single-particle subband 
with respect to fi. Figures [2](a), (b) and (c) show ipnm{0,z) [given in units of V'nm(0,0)] as 
a function of z for {n,m) = (0,0), (0,±1) and (0,±2), respectively (here a = — 140nm). 
Curves for fi/huj_ = 1.8, 2.2, 2.8 and 3.1 are given in each panel. For all these values of 
fi/huj_ the bottom of single-particle subband (0,0) [see panel (a)] is situated far below the 
Fermi surface and, as a result, ipofi{0,z) exhibits features typical for loosely bound Cooper 
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pairs (see, e.g., Ref. 1371) . i.e., we observe an exponentially decaying curve with superimposed 
fast oscillations with period 2TT/kp'^\ where k^p'"^^ = ■\J ^ [/i — huj±{l + 2n + \m\)] is the 
subband-dependent longitudinal Fermi wavevector. Unlike Fig. |2]^a). Figs. |2|^b) and (c) show 
results that significantly differ from the typical BCS behavior [except of fi/hw± = 2.8 and 3.1 
in panel (b), where the bottoms of the subbands with {n, m) = (0, ±1) goes below the Fermi 
surface]. These results resemble the behavior of a wave function of a real bound state in a 
two-body problem unless the bottom of the corresponding subband goes below the Fermi 
level, see fi/hLO± = 1.8 for {n,m) = (0,±1) in panel (b) and fi/hw± = 1.8, 2.2 and 2.8 for 
{n,m) = (0, ±2) in panel (c). We note that for the particle densities in the center of the 
trap Up ^ 2 X 10^^ - 5 x 10^^ cm^'^ (this corresponds to our choice of the chemical potential, 
see details below) the mean inter-particle distance l/n^^ varies in the range 0.6 - 0.8 /im. 
As /|| ~ 2.5/im, than ^g'''^^'' is less than the mean inter-particle distance for fi/hoo± = 1.8 [see 
Fig. El^b)], and the same is true for ^o°'^^'' when jj,/huj± < 2.8 [see Fig. EJ^c)]. Hence, we 
observe a universal behavior (not depending on particular quantum numbers) of ipnmiOjz) 
that exhibits a strong localization when the perpendicular level hjuj±{l + 2n+\m\) crosses the 
Fermi surface. Such a localization of the pair distribution is to a great extent similar to the 
behavior of the fermionic pair-wave function at the ordinary BCS-BEC crossover driven by 
changing the pair interaction strength^!. What is the reason for the longitudinal localization 
of the pair distribution in our case? When the bottom of a single-particle subband is situated 
far below the Fermi level, the ratio of the interaction energy to the longitudinal kinetic 
energy in such a subband is small. As a result, we obtain the ordinary BCS picture for 
the longitudinal distribution of fermions in a condensed fermionic pair. However, when the 
bottom approaches the Fermi level, the longitudinal Fermi motion is depleted, resulting in a 
significant decrease of the ratio of the longitudinal kinetic energy to the interaction energy 
and, as a consequence, in the longitudinal squeezing of a condensed fermionic pair. The 
depletion of the longitudinal Fermi motion in a subband with the bottom close to fi can be 
seen from Figs. Mjo): the period of oscillations in ^/'o,±i(0, -z) decreases when passing from 
/i = 2.2hw± to 3.1^_L. Here the question can arise if a subband with the bottom above the 
Fermi level can contribute to the pair condensate. For an ideal Fermi gas such a subband is 
not populated. However, this is not the case in the presence of superfiuid/superconducting 
correlations that smoothen the Fermi surface. If the energy spacing between the bottom of 
a subband and the Fermi surface is about or less than the relevant pairing gap, this subband 
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still contributes. In the opposite case its role is diminished and becomes negligible with 
an increase of the above spacing. For instance, as seen from Figs. [H^b), (d) and (f) the 
subbands with (n, m) = (0, ±1) make a significant contribution to the coherent properties 
even at /i = 1.9^^. However, their contribution becomes negligible for fi < 1.8hw±, where 
the total longitudinal length approaches ^q'''^''. It is also of importance to note that the 
ratio of the interaction energy to the total kinetic energy is not affected by the perpendicular 
quantization, i.e., it remains small enough even when the bottom of one of the single-particle 
subbands approaches the Fermi surface. The effect of interest is due to a redistribution of 
the kinetic energy between the parallel and perpendicular degrees of freedom in the subband 
whose bottom approaches fi. 

In the previous paragraphs we considered the impact of the multi-subband structure on 
the off-diagonal longitudinal behavior of \l/(r, r'). It is of interest to check what happens with 
the diagonal function \l/(r, r) that controls the order parameter A(r) = 5'\E'(r, r) and, so, the 
spatial distribution of the fermionic condensate. In Fig.[3]the contour plots of A(p, z) (calcu- 
lated for T = and given in units of fi) are shown together with the partial contributions of 
the relevant subbands at a = — 140nm for fi/huj± = 1.9, 2.1 and 2.8. The left contour plot of 
each panel shows the order parameter whereas the two right contour plots depicts the partial 
contributions to A(r) of the subbands with (n, m) = (0,0) [lower] and (n,m) = (0, ±1) [up- 
per]. At fi/fioj± = 1.8, see panel (a), the quantum-size resonance associated with the two 
degenerating subbands (0, ±1) begins to develop. So, the longitudinal distribution of the 
order parameter is mainly determined by the states with (n, m) = (0,0). As a result, the 
distribution of the condensate distribution is rather extended in the z direction, i.e., approx- 
imately from — 11 /|| to 11 /||. Two points of enhancement of the order parameter can be seen 
next to the left and right edges of the longitudinal condensate distribution: at z/l\\ = —9.3 
and 9.3, respectively. These local enhancements are signatures of the BCS character of 
the pairing in the quantum channel (0,0): their 2;-coordinates are solutions of the equation 
Mu\\z^/2 = h^{kp°^y/2M. Unlike {n,m) = (0,0), subbands with (0, ±1) produce a con- 
tribution localized around z = 0, which is in agreement with the quasi-molecule character 
of the fermionic pairing in these subbands. In panel (a) such a contribution does not have 
a significant effect on the longitudinal distribution of A(r). However, when the bottoms of 
subbands (0, ±1) cross fi, the corresponding shape resonance develops and the condensate 
distribution acquires a clear bimodal character, as seen from panel (b). Here the main con- 
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tribution to the pairing correlations comes from the states with (n, m) = (0, ±1), which gives 
rise to a significant enhancement of the order parameter around z = 0, i.e. from z = — 3/|| 
to 3/||. Such a locahzation along the z direction disappears when the perpendicular level 
2hw± goes significantly below fi, as seen from Fig. ^c). Here we arrive at the BCS picture 
of the fermionic pairing in subbands (0,±1), similar to that of {n,m) = (0,0). 

Next we discuss important physical characteristics such as T^/Tp and kF\a\ that are 
directly related to the problem of cooling the Fermi gas—. Typical values of kp\a\ reached 
at the BCS side of the Feshbach resonance are down to ~ 1 whereas the corresponding 



values of T^/Tp fall into the interval 0.05 - 0.2 (see Ref. |22| ). To check if our choice of is 
consistent with these experimental results, we need to estimate kp, the Fermi wavevector 
in the center of the trap. The simplest way to estimate kp is io extract its value from the 



difference /i — i.e., to use kp = y -^(/^ — hu±). Such an estimate tells us that kp\a\ 
varies from 0.56 to 0.7 at a = — 140nm when ^/hLO± increases from 2 to 4. At the same 
time, as follows from Fig. [TJ^a), TjTp [Tp = ffkj,/2M) varies from 0.06 down to 0.03 (see 
also Fig. 2 in Supplementary Information). For a = — 210nm we have kp\a\ = 0.84 and 1.05 
at ^/hu)± = 2 and 4, respectively. Here T^/Tp goes from 0.14 down to 0.07 (Supplementary 
Information, Fig. 2). A more accurate estimate oi kp can be found when calculating the 
position dependent particle density np(p, z) (see the relevant formulas in Supplementary 
Information). Contour plots of np{p,z) [calculated at zero temperature for a = — 140nm 
and given in units of 10^^ cm~^] are shown in Fig. |4] for fi/hw± = 1.7 (a), 2.1 (b), 2.8 (c), 
3.1 (d) and 3.7 (e). In general, from Fig. H] one can extract particle densities close to 



those of the simplified estimate. Indeed, based on kp = y ^(/x — huj±), one finds that the 
mean particle density is about 2 x 10^^ cm~'^ at fi/hLO± = 2 and increases to 4 x lO^^cm"^ 
when hw± approaches 4. From Fig. H] it follows that the mean density in the center of 
the trap is about 4 - 5 x 10^^ cm~^. However, it does not show any systematic increase 
with decreasing u±. As follows from Fig. HI np{p, z) in the center of the trap oscillates 
when changing co±. Such oscillations are due to an interplay of the relevant single-particle 
subbands whose number changes by one each time when fi/hw± approaches an integer. Now, 
taking Up = 4.5 x lO^^cm"'^ in the center of the trap and using a = — 210nm, we obtain 
that kp\a\ = 1.07 and Tc/Tp varies between 0.07 and 0.09 (see Fig. 3, Supplementary 
Information). Thus, both ways of estimating kp\a\ and Tc/Tp give experimentally attainable 
values. We remark that typical densities of ultracold trapped alkali-metal gases are about 
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10^2 . iQis (see Ref. 



221 ) ■ We would also like to note that the quantum-size oscillations 



of the coherence length are much more pronounced for smaller couplings, i.e., for fci?|a| ^ 1. 
However, for such couplings the ratio T^/Tp becomes significantly below its experimentally 
attainable lower bounds 0.05 - 0.2. 

We also note that small quantum-size oscillations in np(p, z) are both qualitatively and 
quantitatively different from those observed in the fermionic condensate and, so, can in no 
way explain the qualitative changes in the pairing correlations discussed above. Although the 
common origin of oscillations in both single- and many-body characteristics is the formation 
of the subband structure, the superfiuid quantum-size oscillations cannot be understood by 
looking at the single-particle density np{p,z). For example, when the system goes through 
the first resonance, i.e., 1.7 < fi/fiuj± < 2.8, the single-particle density in the center of the 
trap changes by about 20%, see Fig. HJ Based on the homogeneous BCS theory, one can 
estimate that the corresponding change in the BCS coherence length is only about 7%, which 
is far smaller than observed in our calculations for a = — 140nm. 

Another interesting quantity to discuss is the product kp^o that is used to to check the 
evolution of the fermionic condensate at the BCS-BEC crossover. According to the paper—, 
the intermediate region of a BCS-BEC crossover is approached when l/vr < kp^Q < 2it. 
This criterion corresponds, for a uniform 3D system, to the domain — 1 < (kpa)^^ < 1, as 



discussed in Ref. 
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The upper boundary can be rewritten as C,o = Xp. When < Xp, a 
condensed fermionic pair has an average size of the order of or less than the interparticle 
spacing. Such a pair acquires a local character: it weakly, or even does not, overlap with other 
fermionic pairs. The lower boundary kp^o = l/^r sets the BEC regime (i.e., for kpC,o < I/tt) 
with the formation of stable point-like molecules. Let us now consider how the product kpC,o 
changes with fi/huj± in a quasi-lD fermionic condensate. The most interesting is the case of 
a = —140 nm. Here kp^o is well above 27i except for the points of the size-dependent drops of 
^0, where this product reaches values between 4 and 6 (see Supplementary Information). It 
means that each time when the bottom of a single-particle subband crosses fi, we approach 



the intermediate region of our BCS-BEC-like crossover driven by quantum-size effects. As 

50 



for the subband-dependent length ^q^"^\ it is always well below 27r in the vicinity of the 



corresponding quantum-size resonance. For instance, kpC,^'"^^^ (taken at fi/huj_ = 2) is about 
2.0 - 2.5 (see Supplementary Information). This is of importance because, e.g., for fi/hw± = 2 
about 70 - 75% of the condensed pairs come from the subbands with {n,m) = (0, ±1), 
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as already mentioned above. Note that for the subband-dependent coherence length it 
looks more natural to check the product k^p'"^^C,Q^"^^ that gives significantly smaller values as 
compared to feF^o"*"^- However, in the vicinity of a quantum-size resonance, k^p"^^ associated 
with the resonant subband significantly underestimates the mean longitudinal single-particle 
momentum due to the neglect of the contribution from pairing correlations. Moreover, 
^(nm) defined when the bottom of a subband is positioned above /i. So, a proper 

generalization of the criterion l/vr < kp^Q < 27r for the subband-dependent BCS-BEC 
crossover can be of interest. 

Several remarks are in order about the ultraviolet regularization used in our calculations. 
As mentioned above, to arrive at an ultraviolet finite approach, we adopt the regularization 
by subtracting l/{2Tnmj) in Eqs. ([I]) (taken with k — )• z/ = {n,m,j}), (jl]) and which 
is a simple and straightforward extension of the regularization for homogeneous superfiuid 
systems to our spatially nonuniform case. In fact, our procedure is a simplified version 
of a comprehensive scheme for spatially nonuniform systems suggested and investigated 



in Ref. 
Ref. 
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Such a simplification is very similar to the ultraviolet regularization used in 
43l . We stress that our basic conclusions are not sensitive to the specific form of the 
ultraviolet regularization. In particular, almost the same results are obtained with a simple 
cut-off procedure, i.e., the term —l/{2Tnmj) is removed from Eqs. ([1]), (jl]) and ([5]) and the 
summation in the relevant expressions runs over the single-particle states with \Tnmj—f^\ < yW, 



see Fig. 4 in Supplementary Information. A similar ultraviolet cut-off is often used in 



apers 
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44. 



dealing with superfluidity in nuclei and uniform atomic condensates, see, e.g., refs. 

Concluding this section, we note that the longitudinal shrinking of condensed fermionic 
pairs driven by quantum-size effects in a cigar-shaped fermionic condensate can be investi- 
gated experimentally via rf-spectroscopy^^. Another possibility is to exploit the presence of 
a coherent superposition of the subband-dependent fermionic condensates. When a shape 
resonance develops, one of them has a molecule-like character and the others are of the BCS 
type. Such fermionic condensates will evolve in a different manner when switching off the 
longitudinal conflnement. Supplemented by measurements of the quantum-size oscillations 
of the basic quantities, i.e., Tc and relevant pairing gaps (the gaps can be probed by the 
rf-spectroscopy), such experiments could give useful information about the quantum-size 
driven reconstruction of the fermionic pairing in a confined condensate. 
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II. QUANTUM SUPERCONDUCTING NANOWIRES 



The effect in question is universal and can also be expected for high-quality supercon- 
ducting nanowires, when disorder is small enough and, so, scattering on imperfections do 
not shadow the formation of single-electron subbands due to the transverse quantization 
of electron motion. For conventional superconducting materials, e.g., Al or Sn, the bulk 
(zero-temperature) excitation gap Abuik varies from^^ ~ 0.1 to ~ 1.0 meV and, so, the inter- 
subband energy spacing 5^2^^ (with d the nanowire diameter) is of the same order or 
larger when < 20 — 40 nm. For these diameters quantum-size effects on the supercon- 
ducting properties become of importance. In particular, recently found width- dependent 
systematic shift of Tc to upper values in superconducting aluminum/tin nanowires^'^'^ 
have been attributed to these quantum-size effects^. We remark that at present the nar- 
rowest fabricated superconducting nanowires are aluminum samples with width down to 8 



-11 nm (see Ref. 



20). 



All the formulas for the superconducting condensate in a quantum wire are similar to 
those of the previous section. However the relevant single-electron wave functions will be 
different. The electron motion is not quantized in the z direction (we deal with periodic 
boundary conditions with unit cell of length L) and, so, for ipi{z,z') = ipi{z — z') we have 
(at T = 0) 

where i stands for the set of quantum numbers controlling the motion perpendicular to 
the nanowire; Aj (chosen real) is the subband pairing gap that does not depend on the 
longitudinal wavevector k (due to periodic boundary conditions in the z-direction); and Ei = 
a/ + , where Eik is the single-electron energy measured from the chemical potential, 
i.e., Eik = €i + 1^ — n with Si the discrete perpendicular single-electron level. As seen from 
Eq. ([7j), there is no regularization term similar to that of Eq. ([5]). But instead, the sum in 
Eq. ([7]) runs only over the single-particle states in the Debye window, i.e., with \eik\ < huo, 
where ud is the Debye frequency appearing in the problem due to the phonon mediated 
attraction between electrons. This is why the lower ki and upper k2 cut-off momenta appear 
in the integral in Eq. ([7]). For k = ki, Sik crosses the lower boundary of the Debye window, 
i.e., ki is a nonnegative solution of Sik = —huj^- This is when the subband bottom is 
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situated below jj, — Hujd- When the subband bottom goes above — Hwd, is set to zero. 
In turn, k2 is a nonnegative solution of Eik = fiWD provided that the subband bottom is 
below the upper boundary of the Debye window. When Si > fi + Hud, we get k^p = and, 
so, such a subband does not contribute to the superconducting characteristics. A major 
contribution to the integral in Eq. ([7]) comes from the states in the vicinity of the minimum 
of \eik\ so that in most interesting cases one can simply put ki = and /c2 — +oo (the 
resulting integral is perfectly convergent). It means that the presence of the cut-off does not 
significantly influence the longitudinal properties of a quasi- ID superconducting condensate 
and, in turn, our conclusions. Based on Eq. ([7]), one can calculate \l'(r,r') from Eq. (jlj), 
where n, m are replaced by i and is no longer the eigenfunction of the 2D isotropic 
harmonic oscillator but the single-electron wave function corresponding to the transverse 
quantum number i. 

When inserting ki = and k2 = +oo into the integral in Eq. ([7]), the decay length of 
'ipnmiz — z') can be calculated analytically by using the contour integration in the complex 
plane (this length is proportional to the longitudinal subband-dependent coherence length 
^Q*^). The integrand in Eq. (I7j) has four singular points (the square-root branch points) with 
the same absolute value for the imaginary part. For, say, positive z — z', the contour should 
be closed in the upper half plane and distorted to encircle the cut between the two upper 
singular points having the same imaginary part [me(A//if~+~^~/^i)]^^^/^) where /ij = ^i — Si 
is the chemical potential measured from the subband bottom. This imaginary part controls 
the decay of ipi{z — z') with increasing z — z' , so that we obtain 

-1/2 

(8) 

When ^i/ ^ 1, i.e., the ratio of the longitudinal kinetic energy to the interaction energy 
in the corresponding subband is large enough, we arrive at the conventional result for the 
BCS coherence length, i.e., ^q*'' oc /ifj/Aj, with Vi = A/2/ii/me the longitudinal subband- 
dependent Fermi velocity. This is for a subband with the bottom far below the Fermi 
surface so that Vi ~ vp, with vp = ^l^ijmf,. As Aj oc Tc, we get oc hvp/Tc-, which 
is in agreement with our results for ^o°'°^ of ^ cigar-shaped fermionic condensate in Fig. [H 
At the point of a quantum-size superconducting resonance, when /Xj — )■ 0, Eq. ^ reduces 
to a completely different expression, i.e., ^'^^ ~ fij {vrif^l^^^l'^ . This is also fully consistent 
with our numerical results given in Fig. [1] (see the discussion in the previous section about 
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n 



xig™' for (n, m) = (1,0), (2,0) and (3,0) at points fi/hu!± = 2, 3 and 4, respectively). When 
fii < and \fii\ ^ Aj, one obtains ^ oc h/{2me\fii\y^'^. Here the dependence of ^g*"* on the 
subband pairing gap Aj (and, so, on Tc) disappears entirely. The same was found in the 
previous section for ^g^'^'' for fi/huj± < 3.5, see Fig. [TJ 

The simple analytical structure of Eq. highlights a strong similarity of our results with 
the BCS-BEC crossover driven by a change of the strength of the pair interaction: /ij and Aj 
in Eq. ([8]) can simply be replaced by the chemical potential fi and the gap A to reproduce 
the evolution of the size of the condensed fermionic pairs through the region of the crossover 
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However, there is also an important 



driven by the Feshbach resonance (see, e.g., Ref. 
difference: in our case a trend similar to a BCS-BEC crossover is found for the resonant 
single-particle subband whose bottom is close to the Fermi surface. Other relevant subbands 
remain in the ordinary BCS regime. Thus, as already mentioned above, at the point of a 
quantum-size resonance the pair condensate is governed by a coherent superposition of the 
quasi-molecule quantum channel with a set of ordinary BCS channels. Resonant channels 
significantly contribute to the condensation energy. For example, this is about 70% and 50% 
for a superfluid cigar-shaped Fermi gas at fi/f}LO± = 2 and 3, respectively (see the previous 
section). For superconducting cylindrical nanowires with diameters less than 8 — 10 nm 
the contribution of the resonant subbands at a quantum-size superconducting resonance 
is typically about 60 — 70%. This makes it possible to conclude that sufficiently narrow 
superconducting nanowires at the resonant points are mainly governed by the quasi-molecule 
channel of the fermionic condensate. As an illustration. Figs. [5] (a) and (b) show contour 
plots of the off-diagonal superconducting order parameter g'^{p, v^, z; p, ip, z') (the transverse 
coordinates of two electrons are taken the same) for a cylindrical superconducting nanowire 
at diameters d = 3.85 nm [panel (a) is for the case of a quantum-size superconducting 
resonance] and d = 4.08 nm [panel (b) is for the non-resonant case]. Such a representation 
is very convenient because it provides us with information of not only the longitudinal 
correlations but, in addition, on the (diagonal) order parameter A(p) = (7\I^(r,r). Unlike 
the previous case of the cigar-shaped trap, the order parameter does not depend on z due 
to periodic boundary conditions in the z direction. As seen from Figs. |5]^a) and (b), the 
longitudinal distribution of electrons in a condensed pair significantly shrinks at d = 3.85 nm 
as compared to d = 4.08 nm. This is because ai d = 3.85 nm the bottoms of the two 
single-electron subbands with the quantum numbers i = {n,m) = (5, ±1) and (1,±11) are 
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situated in the vicinity of the Fermi surface (due the cyhndrical shape the relevant transverse 
quantum numbers are again n and m, the radial and azimuthal numbers, respectively). 
The corresponding quantum-size superconducting resonance is responsible for a significant 
enhancement of the order parameter A(p) as compared to its bulk value (for our parameters 
it is = 0.25 nm, see Supplementary Information) and for the main contribution of these 
subbands to its averaged value (about 70 %). When increasing the diameter, these subbands 
shift down in energy, the corresponding resonance decays and, so, one can see a rather 
extended distribution along the z direction aX d = 4.08 nm, see panel (b). Fig. (I5])(c) 
demonstrates the longitudinal-decay length of A(p, z — z') calculated from our results by 
a numerical exponential fit for p/R = 0.18. As seen, the variation of the length (it is 
proportional to the longitudinal pair size .^o) for diameters from d = 3.85 nm to 4.08 nm 
is two-three orders of magnitude. It is significantly more pronounced than the drops in 
^0 found in the previous section for an ultracold quasi- ID Fermi gas. The reason is as 
follows. From our analytical results it follows that the resonance- driven decrease of the 



longitudinal coherence length is governed by a factor ^ kkp/VmA, where A is the typical 
energy gap; m is the relevant particle mass (e.g., m = rrie for electrons and m = M for 



lithium atoms); and kp can be estimated as y -^(/i — hu±) for a cigar-shaped Fermi gas 
whereas for a metallic nanowire we have kp = ^J^2/m^W. Then, one can find that the 
reason for the significantly more pronounced drop of the longitudinal BCS coherence length 
in the superconducting nanowire is the value of kp which is more than three orders of 
magnitude larger in metals. The effect of an increase of the typical excitation gap in metals 
by seven orders of magnitude as compared to a trapped atomic condensate is strongly 
weakened by the lO'^-drop of the particle mass and, in addition, by the presence of the 
product mA under the square root. For instance, one can simply compare ^g*^ = hvi/Ai 
with ^Q*'* = /i/(meAj)^/^ for typical metallic parameters of weak-coupling superconductors: 
f j ~ 1 X 10^ — 2 X lO^m/s and Aj ^ 0.2 meV (see. e.g., Ref. |45|). This allows one to find 
that ^0*^ drops at a resonant point by two-three orders of magnitude, from the micron to 

(i) 

the nanometer scale, i.e., down to ~ 20 nm. Then, taking account of the fact that the 
diagonal order parameter A(p) = A(p, z — z')\z-z'=o is 3 — 4 times enhanced as compared 
to bulk aX d = 3.85 nm, one can eventually obtain a drop shown in Fig. [5] (c). We note that 
when moving below diameter 3.85 nm, the longitudinal decay-length first increases with 
decreasing d but then, when reaching d = 3.82 nm, it begins to drop. This is the effect of 
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another quantum-size resonance situated at d = 3.72 nm (see Ref. |32|). For diameters larger 
than 3.85 nm, the next quantum-size resonance develops at d = 4.2 nm. 

We would like to remark that the above results for superconducting nanowires were 
calculated in the clean limit. In the presence of surface imperfections and disorder of real 
metallic nanowires the quantum-size oscillations of the longitudinal coherence length will 
be smoothed. However, such giant drops of can hardly be completely washed out (see 
discussion in Supplementary Information). Yet, it will be difficult to observe details of 
oscillations of the BCS coherence length due to fluctuations of the confinement dimensions 
of real samples. To capture the details, it is more promising to use materials with lower 
charge-carrier densities (e.g., semimetals and doped semiconductors). Such materials have 
larger Xp that controls the number of the occupied transverse modes and the width of 
the superconducting resonances (this width can be defined as a maximal variation of the 
nanowire diameter with no significant effect on a resonant enhancement). We can expect 
that for Xp ^ Inm smoothing of quantum-size oscillations due to width fluctuations will 
be less signiflcant. 



III. CONCLUSION 

Since the classical paper by Cooper~ it is well-known that the conflguration of the phase 
space available for scattering of time-reversed fermions plays a crucial role for the formation 
of condensed fermionic pairs. Indeed, only a strong enough attractive interaction between 
fermions with opposite spin in 3D is able to produce a two-body bound state in vacuum. 
However, when scattering of fermions is influenced by the exclusion of a fllled Fermi sea, 
i.e., the available phase space is restricted by exclusion of the single-particle states inside 
the Fermi sea, we arrive at the Cooper instability resulting in the formation of weakly 
bound in-medium pairs of fermions for arbitrary strength of the attractive interaction. Re- 
stricting the phase space by removing the fllled Fermi sea, one actually removes long range 
contributions to the Cooper-pair wave function, which, say, "encourage" fermions to form 
in-medium bound states. Our present results show that the additional reconflguration of the 
phase space due to quantum conflnement, such that the band of single-particle states splits 
up into a series of subbands, can further modify the scenario of pairing. The formation of 
multiple subbands due to quantum conflnement results in a signiflcant redistribution of the 
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kinetic energy between confined and unconfined directions. In particular, for a quasi-lD 
fermionic condensate tlie subband-dependent ratio of the longitudinal kinetic energy to the 
interaction energy drops down to almost zero when the bottom of a subband crosses the 
Fermi level (this is rather similar to superconducting semimetals with the chemical potential 
driven by superconducting correlations below the bottom of the conduction band^). The 
longitudinal Fermi motion of fermions in such a subband is depleted, which results in a drop 
of the longitudinal size of condensed fermionic pairs. So, each time when the bottom of 
a single-particle subband passes through the Fermi level, a quantum-size resonance devel- 
ops and the sup erfluid/ superconducting system exhibits trends similar to the well-known 
BCS-BEC crossover driven through the Feshbach resonance in ultracold Fermi gases in 3D 
traps (see, e.g., Ref. |22i) . In this case condensed fermionic pairs behave, to a great extent, 
similar to condensed boson-molecules. For instance, this is reflected in a clear bimodal spa- 
tial distribution of a harmonically trapped quasi-lD fermionic condensate. A contribution 
governed by the quasi-molecule channel associated with the subband whose bottom is situ- 
ated in the vicinity of the Fermi surface is strongly localized in the center of a trap. Other 
quantum channels are due to single-particle subbands with bottoms far below the Fermi 
surface. They are in the BCS regime and yield an extended longitudinal distribution of the 
fermionic condensate typical of Cooper pairs. When the number of relevant single-particle 
subbands increases, the role of quantum-size effects is diminished and, finally, all related 
phenomena disappear. One can expect that the impact of the subband formation on super- 
fiuid/superconducting properties is of importance when the inter-subband energy spacing is 
larger than or close to ~ fcsTc. We remark that similar physics can be expected for quasi-2D 
fermionic condensates. 
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Supplementary information 

Appendix A: Quasi-ID superfluid Fermi gas 

Our investigation is focused on the quasi- ID Fermi gas, where the single-particle energies 
Tnmj (with j corresponding to the longitudinal motion) with the same n and m form a 
quasi-continuum. So, it is convenient to introduce single-particle subbands (n, m) as seen 
from Fig. [6l In panel (a) we show Tnmj versus the quantum number j for fi/hw± = 2. 
In this case the bottom of subband (0, 0) is located far below the chemical potential 
whereas the bottom of two degenerate subbands (0, ±1) almost touches fi (it is higher only 
by hw\\/2 = O.OOSyU, where we take fi = h ■ 24 kHz, see the discussion about the chosen 
parameters in the article). When decreasing uj_ (and keeping fi the same), more and more 
new subbands crosses the chemical potential. In particular. Fig. [6](b) illustrates the case of 
fi/huj_ = 3. Here the bottoms of subbands (0,0) and (0, ±1) are below /i. The bottoms 
of the three degenerate subbands with (n, m) = (0, ±2) and (1, 0) almost touch /i, which 
results in the appearance of a new quantum- size resonance. 

When discussing in the article important physical parameters of the superfluid fermionic 
condensate such as T^/Tp, fci?|a|, and kp^Q, we invoked two different ways to estimate kp 
in the center of the trap. The simple one involves the chemical potential /i measured 
from the lowest single-particle level hw± + |^|| ~ huj± and is based on the relation 



kp = Y^(/U — ^i^±)- This estimation is good enough when a trap is almost isotropic, 
i.e., the trapping frequencies are close to each other in all directions. For the quasi-lD case, 
when uj± ^ u\\, the reliability of the above estimate can suffer from the anisotropy in the 
distribution of the kinetic energy in single-particle subbands. For instance, when the sub- 
band bottom is situated in the vicinity of /i, the longitudinal kinetic energy is suppressed 
as compared to the kinetic energy of the motion perpendicular to the system. While the 
subband bottom goes far below /i, the above anisotropy disappears in favor of another: now 
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the longitudinal motion prevails. So, the resulting values of kp and Tp = fi^k'^/2me will be 
dependent on whether kp is estimated from the longitudinal kinetic energy or it is extracted 
from the energy of the perpendicular motion. Another issue of uncertainties is the choice of 
the relevant subband to extract kp. Notice that the above simple estimate is based on the 
longitudinal kinetic energy of the lowest single-electron subband. 

A more accurate (but more time consuming) procedure to estimate the relevant kp m 
the center of the trap is to calculate the spatial particle distribution ne{v) and, then, to find 
the mean particle density in the center of the trap. In the mean-field approximation the 
quantity ne(r) is calculated through the coherence factors Uj,{v) and f,y(r), see Eq. (3) in the 
article, i.e., 

n,(r) = 2 (1 - /.) + M)? /.] , (Al) 

V 

with ]y = l/le^^" + 1), the Fermi distribution of the bogolons. We note that for a trap 
with axial symmetry the density of particles depends only on the radial coordinate p and 
the longitudinal coordinate z. Based on Eq. (3) of the article, one can rewrite Eq. ( lAll) as 
(at T = 0) 

'■'nmj 

nm] ■' 

where \'dnm{,P-, V^)! does not depend on the azimuthal angle 99. We note that no regularization 
similar to subtracting the term XjlTnmj in Eq- (5) of the article is now needed. The sum 
in Eq. flA2p is convergent and, in addition, the regularization used in Eq. (5) of the article 
assumes the explicit presence of the coupling constant, which is not the case for the position 
dependent particle density np(r), see Eq. ( lAll) . A numerical solution of the self-consistency 
equation (1) of the article provides us with a set of Aj^'s that should be inserted into the 
above relation for np(p, z). Our numerical results for np(p, 2;) are shown in Fig. 4 of the article 
as contour plots. When checking the mean particle density in the center of the trap, one 
can find that it varies from 4 x 10^^ to 5 x 10^^ cm^^. So, when taking = 4.5 x 10^^ cm~^, 
we can find kp = (37r^np)^/'^ = 5.1 fim~^ and \p = 27^ /kp = 1.23 yum. This can be compared 



n,ip,z) = 5^(1 - Z^^^) |^_(p,^)|2;,|(.), (A2) 



with kp = ^ "^{p — hu±) = 3.99 pm^^ and 4.88 pm~^ at p/hu± = 2 and at 4, respectively. 
Both estimates produce close values for kp in the center of the trap but we obtain different 
trends for the mean density in the center of the trap. According to the simple estimate, Up 
plainly increases with p/fioj±. Whereas Up found from the position dependent distribution 
of the particles slightly oscillates around 4.5 x 10^^ cm~^. 
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Having at our disposal possible values of kp and using our numerical results given in Fig. 
1 of the article, we can find Tc/Tp and kp^o as shown below in Fig. [7] for the simple estimate 
and Fig. [H]for the estimate based on ne(p, z). 

Concluding this section, we would like to note that our results are not sensitive to the 
regularization procedure used in Eq. (1) of the article (taken with k — {n,m,j}) in order 
to remedy the divergent sum (the ultraviolet divergence). For illustration. Fig. H] shows 
numerical results for the critical temperature found when the regularization term l/[2Tnmj) 
is abandoned in favor of a simple cut-off \Tnmj — /^l < n^n^ Anmj), i-e., the sum in Eq. (1) 
is only limited to the single-particle states with energies satisfying this inequality. Note t 
similar cut-off is often used in papers on ultracold atomic gases, see, e.g., Refs. 
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44. 



Though Tc calculated with the cut-off regularization are slightly different as compared to 
the results given in Fig. 1 of the article, this difference practically disappears when making 
a small shift in the scattering length, i.e., a — )■ a + 20 nm, see Fig. [91 Thus, a particular 
way of ultraviolet regularization does not influence our results on quantum-size oscillations 
of Tc, and the same holds for our conclusions about the size- dependent oscillations in the 
BCS coherence length. 



Appendix B: Superconducting metallic nanowires: parameters used in the calcula- 
tions 

Figure |5] of the article shows our numerical results calculated for a cylindrical super- 
conducting aluminum nanowire with the unit cell (periodic boundary conditions) in the z 
direction = 50 fim. The coupling constant g is chosen so that gN{0) = 0.18, with A^(0) 
the bulk density of states at the Fermi level per unit volume and per spin projection (this 
value of gN{0) is typical of aluminum, see, e.g., the textbook—). Effects resulting in the 
width dependent coupling are beyond the scope of the present investigation (such effects do 
not influence our main conclusions). The Debye energy is taken the same as in bulk—, i.e., 
^n/kB = 375 K {huJo = 32.3 meV). The electron band mass rrie is set to the free-electron 
mass in the calculations. The Fermi level is put to 0.9 eV. We note that this is an effec- 
tive Fermi level which is used together with the parabolic band approximation to correctly 
reproduce the period of the quantum-size oscillations. It is a parameter that depends on 
the interplay between the main crystalline directions and the direction of quantum confine- 
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ment (for a more detail, see discussion in Ref. |32| ). The above value of the effective Fermi 
level is justified by a good agreement of theoretical and experimental results for the critical 
temperature in aluminum nanowires^^. However, this particular choice is not crucial for our 
conclusions. We also take into consideration a systematic shift up of the Fermi level with 
decreasing d for d < 5 — 6 nm. Such a shift appears if one keeps the same value of the mean 
electron density when changing d. 



1. Surface imperfections 

Our consideration of the superconducting nanowires is in the clean limit and assumes 
mirror reflections from the boundaries. So the question arises to what extent imperfec- 
tions of real metallic nanowires can influence our results. We focus on recently fabricated 
high-quality superconducting nanowires, where the superconducting state survives down to 
diameters of about 8 — 10 nm without any signatures of suppression of the critical temper- 
ature driven by disorder [see refs.— >2S about Al nanowires made of strongly coupled grains 
and paper— about Sn single-crystalline nanowires). Unlike strongly disordered nanosized 
superconductors, such high-quality superconducting nanowires exhibit a systematic shift-up 
of the critical temperature with a reduction in their cross-section (50% in Al specimens and 
10 — 20% in Sn nanowires). As recently shown^, it is the quantization of the transverse elec- 
tron spectrum that shifts T^. up in such nanowires through quantum-size superconducting 
resonances. Thus, disorder is relatively minor in high-quality superconducting nanowires 
and the effects of the formation of well distinguished single-electron subbands become of 
importance. 

To go in a more detail about a role of imperfections in quantum superconducting 
nanowires, we remark that the surface roughness is here a major disorder mechanism. Indeed, 
in most pa.peTSF^>^ the mean free path is estimated to approximately follow the nanowire 
width i d, i.e., elastic scattering on the boundary imperfections controls the electron mean 
free path. So, we need to clarify whether or not the size-dependent drops of the longitudinal 
BCS coherence length found in the clean limit will be significantly influenced by the surface 
roughness. The first thing to do with this is to compare the Cooper-pair size with the 
electron mean free path. Beyond a resonant point increases up to microns and, so, we 
are in the dirty limit, i.e., ^ ^- However, when approaching a point of the quantum-size 
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superconducting resonance, drops down to values close to the nanowire diameter, see, 
e.g.. Fig. [5] in the article (this is true for nanowires with d < 10 — 20 nm, for larger diameters 
the effect is washed out). As the mean free electron path I is about d, then we obtain 
^0 ~ ^- So, the effect of disorder on the superconducting correlations is not significant: the 
suppression factor in Green's functions is e~^/2^ ~ 1 for r ~ ^0- Thus, the elastic scattering 
on the surface imperfections will not significantly alter our conclusions: at points of the 
size-dependent drops we are close to the clean limit. 

The above reasoning assumes that the density of single-electron states at the Fermi level 
does not change significantly due to surface imperfections. So, the next step is to clarify if the 
single-electron spectrum is seriously influenced by surface imperfections. A size-dependent 
drop of the longitudinal fermionic-pair size occurs due to a subband (subbands) whose bot- 
tom (bottoms) is situated in the vicinity of the Fermi surface (see the article). For such a 
subband the longitudinal motion of electrons is suppressed and the most contribution to the 
single-electron energy comes from the transverse spectrum of electrons. This is exactly the 
reason for the quantum-size driven BCS-BEC crossover. Another consequence of this redis- 
tribution of the kinetic energy between the longitudinal and transverse degrees of freedom 
is long longitudinal wavelengths of the electrons in a resonant subband. The longitudinal 
electron energy in a subband whose bottom is situated in the Debye window is about or 
smaller than the Debye energy. This makes it possible to find that the typical longitudinal 
wavenumber is smaller than k < ^/^vrififJcil^- So, this simple estimations suggests that 
the typical longitudinal wavelength of electrons in a resonant subband is larger than 10 nm. 
This is significantly larger than the characteristic size of the surface imperfections usually 
estimated as being of the order of dimensions of crystalline unit cell (smaller than Inm). 
Hence, the longitudinal part of the electron spectrum in a resonant subband is stable against 
surface imperfections. This is not the case for subbands with bottoms positioned far below 
the Fermi level. Here the electron spectrum can be rather sensitive to the surface roughness 
and the corresponding density of states can even be suppressed at the Fermi level in the 
presence of strong disorder. However, the contribution of such single-electron subbands at 
a resonant point is only of secondary importance: the major effect comes from the resonant 
subbands. 

In addition to the influence of the surface roughness on the longitudinal electron motion, 
one should also keep in mind an uncertainty in the transverse energy due to fluctuations in 
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the nanowire diameter. This uncertainty is able to affect the formation of superconducting 
resonances when it approaches the Debye energy that controls the selection of the single- 
electron states making a contribution to the basic superconducting characteristics. The 
energy spacing between single-electron subbands can be estimated as Let us take 

d = do + Ad, where Ad represents a fluctuation part with values from to 1 nm (this a 
typical scale for the diameter fluctuations in nanowires with width less than 10-20 nm). The 
corresponding fluctuation of the inter-subband energy spacing is estimated as ^^Ad, which 
can be considered as an uncertainty in the transverse electron energy. For the narrowest 
high-quality nanowires with d = 8 — 10 nm, such an uncertainty is about 1 meV, which is 
much smaller than the Debye energy for aluminum /iw/) ~ 32 meV. Even for extremely small 
diameters about 4 nm chosen for a simple illustration in the article (to avoid a discussion of 
many single-particle subbands responsible for the formation of superconducting resonance at 
nanowire diameters ~ 10 nm), this uncertainty is close to lOmeV (when taking Ad = Inm). 
For aluminum this is still three times smaller than the Debye energy. 

Another reason for broadening of the single-electron levels is the hybridization with elec- 
trons of a semiconductor substrate, which is strongly dependent on the fabrication condi- 
tions. However, such a hybridization can be expected to be of importance for specimens 
with width down to a few monolayers like in ultrathin superconducting nanofilms. 

Thus, based on the above discussion we can conclude that our results are quite stable to 
imperfections of real superconducting nanowires. These imperfections will definitely smooth 
quantum-size oscillations of but will hardly avoid the effect of interest. 

2. Fluctuations 

We work in the mean-field approximation and, so, one more point to discuss is the 
effect of fluctuations. It is well-known that fluctuations generally play an important role in 
low-dimensional systems. In superconducting nanowires the main focus is usually on the 
phase fluctuations of the pair condensate: thermally-activated and quantum phase slips, 
see, e.g., the paper— i^. These fluctuations lead to a residual resistance remaining below 
Tc in narrow nanowires, corrupting the superconducting state. Below we argue that the 
BCS-BEC crossover driven by quantum confinement is still pronounced above the nanowire 
diameters, where phase fluctuations proliferate. 
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Effect of thermal fluctuations is usually estimated through the Ginzburg-Levanyuk pa- 
rameter calculated from the conventional Ginzburg-Landau functional. However, for quan- 
tum nanowires of width 5 — 20 nm the conventional Ginzburg-Landau formalism can not be 
applied due to the breakdown of the translational invariance in the direction perpendicular 
to the nanowire. A consequence of such a breakdown is that the order parameter strongly 
varies in the direction perpendicular to the nanowire (see, e.g., Ref. |16!) . The characteris- 
tic length for its spatial variations in this direction is about the nanowire width d, which 
prevents us from using the conventional Ginzburg-Landau formalism which assumes that 
the order parameter varies slowly on the scale governed by the zero-temperature coherence 
length. In addition, the ID Ginzburg-Levanyuk parameter is not relevant either because a 
superconducting quantum nanowire is an essentially 3D system with multiple single-electron 
subbands (more than 10 — 20 subbands even for d = 4nm). Such a multichannel system can 
not be reduced to effectively one-dimensional case, and this is also seen from the fact that the 
order parameter has a nontrivial transverse profile that changes significantly with a change 
of d. To overcome this problem, one should abandon the conventional Ginzburg-Landau 
formalism in favor of its multichannel version. However, a significant complication here is to 
accurately incorporate issues related to the quantum-size driven BCS-BEC crossover (e.g., 
fluctuating Cooper pairs with a nonzero center-of-mass momentum). A simpler option is to 
rely upon available experimental results. They suggest that the Ginzburg-Levanyuk param- 
eter is about 0.1 — 0.15 in superconducting nanowires with diameters ~ 10 nm. In particular. 
Fig. 1 of the paper— (for zero magnetic field) demonstrates that the nanowire resistance falls 
by two orders of magnitude when temperature reduces from T = Tc to T = 0.85 Tc. This 
is still a rather sharp transition with the thermal broadening of about ST/Tf. = 0.1 — 0.15 
and, so, the mean-field treatment looks quite reasonable. 

From the same Fig. 1 in the paper— one can learn that quantum phase fluctuations 
in an aluminum nanowire with width 10 nm result in a residual resistance even below 
T/Tc = 0.5 — 0.6. At this nanowire width such a resistance is smaller than 10~^ in units 
of the normal resistance. It is however expected that for diameters < 10 nm quantum- 
phase slips proliferate, which leads to a superconductor-to-normal crossover at d = dc, with 
dc < 10 nm (see, e.g., papers^^i^). Recent results of the paper— suggest that d^. ^ 8nm. 
Yet, it is rather difficult to analyze experimental data for very narrow nanowires because it is 
not possible to completely rule out weak links as the sources of the residual resistance. Thus, 
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we can conclude that the best regime to probe any signatures of the BCS-BEC crossover 
induced by quantum confinement is to investigate the superconducting nanowires with di- 
ameters just above d^- We would ,like to note that the effect of interest is pronounced for 
d < 10 nm and completely washed out only when d > 20 nm. In addition, the subject of our 
investigation is the zero-temperature coherence length tha t is on e of important parameters 
controlling the rate of quantum-phase slips, see, e.g., Ref. 
expect that our results can be of relevance even for d < dc- 
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This makes it possible to 
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Figure 1 



FIG. 1: Superfluid Fermi gas confined in a cigar-shaped harmonic trap: the critical temperature 
Tc (in units of ^) and the BCS coherence length ,^0 (in units of /||) versus ^/hw^. Panels (a) and (b) 
for a = — 140nm; (c),(d) for a = — 180nm; (e) and (f) for a = — 210nm. The subband-dependent 
coherence lengths ^(""^^ for (n,m) = (0,0), (0, ±1) (0, ±2) and (0, ±3) are also given to compare 
with ^0- 
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Figure 2 

FIG. 2: Cigar-shaped superfluid Fermi gas: the absolute value of the subband-dependent longitu- 
dinal fermionic pair "wave function" (0, 2:) /V'(0, 0) (i.e., one fermion is fixed at the origin and 
another is at the point z) versus z/l^^ for a = — 140nm at fi/huj± = 1.8, 2.2, 2.8 and 3.1. Panels 
(a), (b) and (c) show the results for {n,m) = (0,0), (0,ibl) and (0,ib2), respectively. 
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Figure 3 

FIG. 3: (a) Contour plots of the order parameter A(p, z) [in units of ^, the left panel] together 
with the partial contributions to the order parameter of the subbands with (n, m) = (0, 0) [the 
lower right panel] and (0, ±1) [the upper right panel] for a cigar-shaped superfluid Fermi gas at 
li/hu}± = 1.9 at a = — 140nm; (b) and (c) show the same but for n/huj± = 2.1 and 2.8. 
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(a) /V/!<w^ = 1.7, a = -140 nm 




Figure 4 

FIG. 4: Contour plots of the position dependent particle density np{p,z) (in units of 10^^ cm^^) 

for fj,/huj± = 1.7 (a), fi/Huj± = 2.1 (b), n/h± = 2.8 (c), fj,/huj± = 3.1 (d) and fi/huj± = 3.7 (e). For 

all the panels the scattering length is a = — 140nm. 
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Figure 5 



FIG. 5: Superconducting nanocylinder: (a) contour plot of g^{p,ip, z; p,ip, z') (given in units of 
the bulk gap = 0.25 meV) in the presence of the superconducting resonance at d = 3.85 nm; 
(b) the same as in the previous panel but now for the non-resonant diameter d = 4.08 nm; (c) the 
longitudinal decay length for the distribution ^'(/9, if, z; p, cp, z') (it is proportional to ^o) versus the 
diameter near the resonance point d = 3.85 nm. 
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FIG. 6: The single-particle energy T^^j versus the longitudinal quantum number j for different 
subbands {n,m) at fj,/tlu)_i_ = 2 (a) and fj./huj± = 3. 
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FIG. 7: Tc/Tp and A;i?^o as calculated from the results given in Fig. 1 of the article on the basis 
of the simple estimate of kp = \J'^{n — hoj±) (with Tp = h^k'^p/2M). Panels (a), (b) represent 
a = — 140nm; (c) and (d) are for a = — 180nm; (e), (f) show the data for a = — 210nm. The 
quantities kp^^^"^^ for {n,m) = (0,0), (0,±1) (0, ±2) and (0, ±3) are also plotted to compare with 
kp^Q. 
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FIG. 8: The same as in the previous figure but now for kp extracted from the mean particle density 
in the center of the trap, i.e., Up = 4.5 x 10"*^^ cm^^ {kp = (Svr^np)^/'^). 
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FIG. 9: The critical temperature Tc (in units of ^u) versus fj,/huj± for a cigar-shaped superfluid 
Fermi gas of ^Li: circles are results for a = — 120nm found from the BCS-like self-consistency 
equation with the simple cut-off regularization \Tnmj — A*| < A*; squares represent the previous data 
from Fig. 1(a) of the article. 
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